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ABSTRACT 


The  effect  of  friction  on  the  propagation  of  long  internal 
waves  in  a rotating  ocean  is  investigated  theoretically  for  the  simple 
case  of  a two-layer  system  subjected  to  a constant  eddy  viscosity. 
Waves  of  inertial  period  are  excluded.  Typical  distances  for  a 5055 
decrease  in  amplitude  are  found  to  be  1,000  to  2,000  km. 
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Cartesian  coordinates  of  position,  with  z < 
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t: 
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Value  of  j at  time  t = 0 and  position  x « i 
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Acceleration  due  to  gravity. 
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Density  difference  between  two  layers. 
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INTRODUCTION 


Evaluation  of  oceanographic  data  is  at  best  uncertain  where  the 
influence  of  long  internal  waves  cannot  be  determined.  This  is  clearly- 
shown  by  the  California  current  where  an  apparent  eddy  structure  is 
found  which  is  caused  by  internal  oscillations  of  diurnal  period 
(Defant,  1950a),  From  the  apparent  wide  distribution  of  such  internal 
waves  it  is  reasonable  to  suppose  that  the  possibility  of  error  exists 
in  the  interpretation  of  much  similar  data,  A knowledge  of  the  occurrence 
and  behavior  of  long  internal  waves  in  the  ocean  becomes  then  of  prime 
importance  for  the  advance  of  oceanography. 

Obtaining  information  on  internal  waves  at  sea  is  a difficult 
process  which  requires  long  series  of  observations  for  statistically 
valid  results.  It  has  been  demonstrated  that  reported  internal  waves, 
in  many  cases,  may  be  only  harmonic  representations  of  short  series  of 
random  data  (Haurwitz,  1952,  1953}  Rudnick  and  Cochrane,  1951),  In 
addition,  in  order  to  adequately  describe  the  behavior  of  the  internal 
waves  it  is  necessary  to  take  measurements  simultaneously  at  three 
stations  (Haurwitz,  1953), 

The  essential  reason  for  the  difficulty  in  obtaining  internal 
wave  data  in  the  ocean  is  the  lack  of  an  adequate  theory  to  guide  the 
investigations.  Even  the  presence  of  tidal  internal  waves  is  not  well 
explained*  The  resonance  theory  (Haurwitz,  1950j  Defant,  1950)  is 
inadequate  because  of  the  large  frictional  damping  and  narrow  belt  of 


resonance  for  these  waves.  Since  internal  waves  of  tidal  period  have 
not  generally  been  found  in  the  open  ocean  but  mainly  in  offshore 
regions,  a mechanism  of  generation  associated  with  the  ocean  boundaries 
is  suggested  (Rudnick  and  Cochrane,  1951).  Model  experiments  have 
shown  the  generation  of  internal  waves  by  the  passage  of  surface  waves 
over  bottom  irregularities  (Zeilon,  19310.  At  that  time  it  was  felt 
this  did  not  offer  an  explanation  for  the  observed  tidal  internal  waves 
since  their  amplitudes  would  be  decreased  rapidly  by  friction  and 
divergence.  However,  internal  waves  generated  by  the  action  of  surface 
tides  over  the  continental  shelves  will  not  be  subject  to  radial 
divergence  as  from  a point  source,  but  will  undergo  convergence  on 
traveling  seaward  from  a boundary  source.  It  is  proposed  therefore  to 
derive  expressions  for  the  effect  of  friction  on  internal  wave  motion 
for  the  case  of  a rotating  earth. 

MATHEMATICAL  FORMULATION 

For  mathematical  simplicity  the  problem  considered  is  the 
propagation  of  long  internal  waves  through  a two-layer  fluid  on  a 
rotating  disk  subject  to  constant  eddy  viscosity.  A right-handed 
coordinate  system  is  used  with  the  origin  at  the  mean  level  of  the 
interface,  the  z-axis  vertical,  positive  upwards,  and  the  x-axis 
positive  in  the  direction  of  wave  travel.  A primed  quantity  refers 
to  the  upper  layer  and  a double-primed  quantity  to  the  lower  layer. 

The  following  assumptions  are  madet  the  vertical  accelerations  are 
negligible  permitting  use  of  a hydrostatic  pressure  distribution; 


-2- 


all  quantities  are  independent  of  position  in  the  y direction ; the  water 
is  incompressible;  and  u,  v,  w and  variations  in  p and  p are  sufficiently 
small  that  terms  of  higher  order  can  be  neglected* 

The  equations  of  motion  then  become,  for  the  upper  layer: 
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and  for  the  lower  layer: 


ia!L  - 2 (OvM  - V 
<H  ^)z2 


- g 


ip.  _M” 

p 


a* 


(3) 


111.  + 2(0  u"  - v -J& 

dt  dz* 


0 . 


(h) 


The  equation  of  continuity  is 
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for  the  lower  layer. 

The  boundary  conditions  arej 
at  the  free  surface  there  is  no  tangential  stress,  or 


at  2 “ hf, 

au'  . 

4v* 

4z 
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at  the  bottom  there  is  no  slippage,  or 


at  z ■ -h",  u"  « v"  - 0, 


(9) 


Equations  (1)  through  (9)  are  solved  in  a manner  similar  to  that  applied 
to  tidal  waves  on  the  Siberian  shelf  by  Sverdrup  (Sverdrup,  1926), 


SOLUTION 


A solution  is  found  by  assuming 


^ VYxei(at“kx)  , 


" -Yx  i(at-kx) 
c,Q  e e , 


where  Zqi  Y,  o,  and  k are  all  constants  characteristic  of  the 

motion.  The  general  solutions  of  equations  (1)  through  (k)  are  then 
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pA  and  P2  a1"6  constants  given  by  j~\  ~ ~ ana  j ~ “ > 

respectively.  Relations  for  the  constant  coefficients  A,  determined 
from  the  boundary  conditions,  after  some  rearranging  are: 
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where  e ■ -SB-.  a small  quantity  for  internal  waves.  Solution  of 
these  equations  yields  the  following  expressions  for  the  coefficients: 
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The  equations  of  continuity  can  now  be  written 
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Integration  of  the  equations  of  continuity  and  substitution  of  the 
values  for  the  constants  A,  yields  the  following  two  expressions  in 
( Y + ik)  and  e : 
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These  equations  in  their  present  form  are  somewhat  intractable.  In  order 
to  determine  the  relative  importance  of  the  various  terms,  the  dimension- 
less quantities 
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are  introduced.  Equations  3h  and  35  can  now  be  written  in  the  form 
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These  equations  can  be  solved  for  any  given  values  of  a,  co?  v,  h1, 
h",  and  Ap  /p  • However,  a particular  numerical  solution  will  not  indicate 
the  dependence  of  k,  Y , K,  ^ , and  %"  on  the  above  independent  variables. 
For  the  most  important  range  of  variables  it  is  possible  to  approximate 
the  relations  to  give  functional  relationships  between  the  above  quantities. 

The  following  ranges  of  variables  will  be  considered,  using 
c.g.s.  units: 

h’  >7  5 x 103  , 

5 x io* n,  h"  yy  10*  , 

1 ^ V 103  , 

^4  io-3, 

0)^  .725  x 10_li  , 

0.7  x 10“U<  a ^ 3x  10“^  , 
with  the  additional  restriction  that 

h»  y 1.2  . (Uoa) 
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The  derived  constants  then  have  the  following  ranges  * 
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Under  these  conditions  the  following  approximations  are  valid  (to  within 
10  per  cent)* 
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The  unknown  quantities  are  now  split  into  real  and  imaginary  parts s 
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the  unknowns  when  v-»0.  Upon  substitution  of  these  expressions  into 

equations  (hi)  and  (U2 ) and  replacement  of  -2-  by  P,  JL.  by  Q,  the 
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These  equations  can  be  solved  by  approximation  in  a straightforward  manner 
if  Sj_j  Tj_,  and  ^ are  all  sufficiently  small  so  that  powers  higher 
than  the  first  can  be  neglected,  to  give 
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Now  the  following  expressions  can  be  written  for  the  unknowns* 
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In  most  cases  of  interest  the  terms  in  equations  (U7)  and  (i;8)  containing 

2 

■AP  will  be  negligibly  small  and  jli  . . Equations  (U7) 

and  (U8 ) then  can  be  written* 


2 

(Si  + iTx)  - (P  + Q)  (1  - i)  , (51) 

(#1  + i?l)  " P <p  “ Q)^  - i)  • (52) 


Approximate  solutions  of  these  equations  for  the  wave  number,  K,  and  the 
damping  coefficient,  Y,  of  the  internal  waves  ares 


i b2  - ko2)(h>+  h") 

(h»+  h")  vl/2  [ (a  - 2w)3^2  + (a+  2o)3^l 

1 g h*  h"  Ap/p 

[,  U \f2  h«h«  a(o2  - Uo2)1^2 

'1/2  ku  h» 


2gAp/p/  h*h» 


(53) 


. 3/2  3/  2 

( CT  - 2(0  ) + (ct  + 2(j  ) 

lia 


m 


The  solution  for  e is  then* 


- iP. 


h" 


P h»+  h" 


hi 

h" 


1 - ^ (l-i) 


(2V 

i4.a(h*  + h")  (ap  - Uu?)V2 


(55) 
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DISCUSSION  OF  RESULTS 


The  wave  number,  k,  is  increased,  thus  wave  length  and  velocity 
are  decreased  by  the  effect  of  eddy  viscosity.  The  correction  term  in 
the  expression  for  the  wave  number  varies  as  the  square  root  of  the  eddy 
viscosity  and  inversely  as  the  depth  of  the  shallower  layer.  Its  mag- 
nitude is  small  unless  the  restriction  of  equation  (IjOa)  is  violated 
by  approach  to  the  inertial  period. 

The  damping  coefficient,  y t varies  directly  as  the  square  root 
of  the  eddy  viscosity,  and  inversely  as  both  the  square  root  of  the 
density  difference  and  the  three-halves  power  of  the  depth  of  the 
shallower  layer.  It  increases  approximately  as  the  square  root  of  the 
frequency  but  has  only  small  dependence  on  latitude  in  the  range  under 
consideration. 

The  ratio  of  free  surface  to  interface  displacement  is  decreased 
in  magnitude  by  friction  with  a decrease  in  time  lag  for  the  free  surface 
motion.  Again  the  correction  terms  vary  directly  as  the  square  root 
of  the  eddy  viscosity  and  are  small  unless  the  inertial  period  is 
approached. 

Representative  values  for  k and  y under  various  typical  con- 
ditions are  given  in  Table  1. 

The  distribution  of  internal  waves  will  depend  very  strongly 
on  the  value  of  the  damping  coefficient  y . This  quantity  varies  from 
a maximum  of  6,0  x 10”^  to  a minimum  of  0.0li5  x 10”^  for  the  conditions 
considered.  That  is,  wave  amplitudes  would  be  decreased  by  $ 0 % in  a 
distance  of  120  km.  and  15,000  km,,  respectively.  The  actual  distance 
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of  wave  travel  will  depend  strongly  on  the  prevailing  oceanographic 
conditions  and  the  wave  period  but  not  on  the  latitude  within  the  restric- 
tions of  equation  (UOa). 

Let  us  consider  a typical  example.  Since  the  depth  of  maximum 
stability  is  that  of  maximum  velocity  shear,  the  effective  eddy  viscosity 
is  taken  to  be  1 cm.^sec"1.  An  average  value  of  100  m.  is  used  for  the 
upper  layer  depth  and  10”^  for  the  density  difference  across 

the  interface.  Then  at  latitude  l£°  N,  k and  y will  be,  respectively, 

,15  km-1,  and  O.S>0  x 10”^  km”'*’.  for  a 12  hr,  period  and  ,05>5>  km--*-,  and 
0.32  x 10“3  km~\  for  a 2\±  hr.  period.  The  12  hr.  and  2)\  hr.  waves 
have,  respectively,  wave  lengths  of  1*2  km.  and  110  km.  and  distances 
for  a $ 0 % amplitude  decrease  of  l,b00  km.  and  2,100  km.  Figure  1 
shows  the  variation  of  y with  the  depths  of  the  upper  and  lower  layers. 
These  results  are  the  right  order  of  magnitude  to  give  regions  of 
strong  internal  waves  surrounding  bottom  irregularities  and  yet  not 
large  enough  to  give  standing  waves  over  the  entire  ocean. 


DAMPING  COEFFICIENT, 


TABLE  1 


Values  of  wave  number  and  damping  coefficient  under  varying  conditions. 

i 

I 


L 

(degrees ) 

T 

(hrs) 

h» 

(m) 

h" 

(m) 

V 

(cm^sec"^) 

Ap  1 
(gms.cnT-3) 

k -i 

(km”1) 

(knT1) 

50 

12 

100 

1000 

100 

10‘3 

.10 

6.0  x 10"3 

30 

12 

100 

1000 

100 

10“3 

ah 

5.U  x 10“3 

15 

12 

100 

1000 

100 

io“3 

a5 

5. 0 x 10-3 

15 

21 

100 

1000 

100 

10-3 

.058 

3.2  x 10-3 

15 

2b 

100 

1000 

10 

10-3 

.056 

1.0  x 10“3 

15 

2h 

100 

1000 

1 

10“3 

.055 

0.32  x 10“3 

15 

2k 

5oo 

1000 

1 

10“3 

.027 

0.0h5  x 10~3 

15 

2k 

100 

1000 

1 

2 x 10“3 

.039 

0.23  x 10-3 

0 

12 

100 

1000 

1 

10-3 

.15 

0.h9  x IQ”3 

0 

2k 

100 

1000 

1 

IQ’3 

.076 

0.35  X 10“3 
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